An exact solution for the KPZ equation with flat initial conditions 
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We provide the first exact calculation of the height distribution at arbitrary time t of the con- 
tinuum KPZ growth equation in one dimension with flat initial conditions. We use the mapping 
onto a directed polymer (DP) with one end fixed, one free, and the Bethe Ansatz for the replicated 
attractive boson model. We obtain the generating function of the moments of the DP partition sum 
as a Fredholm Pfaffian. Our formula, valid for all times, exhibits convergence of the free energy (i.e. 
KPZ height) distribution to the GOE Tracy Widom distribution at large time. 
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The continuum Kardar-Parisi-Zhang (KPZ) equation 
is the simplest equation describing the non-equilibrium 
growth in time t of an interface of height h(x, t) in pres- 
ence of noise [H. In one dimension, i.e. x € R, it reads 



d t h = isV 2 h+ -A (V/i) 2 



i](x,t) 



(1) 



where rj(x,t)r)(x,t') = DS(x — x')5(t — t') is a centered 
Gaussian white noise. Originally conceived to describe 
growth by random deposition and diffusion, it defines a 
universality class believed to encompass an astounding 
variety of models and physical systems Q. The growing 
KPZ interface becomes, at large t, statistically self-affine 
with universal scaling exponents. In d = 1, its width is 
predicted j3| to grow as Sh ~ t 1 ^ 3 , as observed in ex- 
periments [4, The KPZ problem also maps to forced 
Burgers turbulence Q, and to the equilibrium statisti- 
cal mechanics of a directed polymer (DP) in a random 
potential, the simplest example of a glass [7[ with applica- 
tions to vortex lines Q , domain walls @ , and biophysics 
10]. However, despite its importance and universality, 
the KPZ equation has vigorously resisted analytical so- 
lutions. 

Progress in analytical understanding of the KPZ class 
in d = 1 came from exact solutions of a lattice DP model 
at zero temperature II ill , discrete growth models such as 
the PNG model 12j, 1 131. a symmetric exclusion models 
[T3 | and vicious walkers [HI . An analogous to the height 
field h(x,t) was identified and, in the large size limit, its 
one-point (scaled) probability distribution was shown to 
equal the (scaled) distribution of the smallest eigenvalue 
of a random matrix drawn from the famous Gaussian en- 
sembles, the so-called Tracy Widom (TW) distribution 
[po| . which appears in many other contexts [13] ■ It was 
found [HI, EH that one gets either the TW distribution 
F2(s) of the Gaussian unitary ensemble (GUE) or Fi(s) 
of the Gaussian orthogonal ensemble (GOE) for droplet 
and flat initial conditions respectively. The correspond- 
ing many point distributions were identified as determi- 
nantal space-time processes, the Airy process, Aii for 
droplet, and Ai\ for flat naturally expressed with the use 



of Fredholm determinants [19( . 

These advances gave valuable, but only indirect infor- 
mation on the continuum KPZ equation, i.e. a conjecture 
for its infinite t limit (termed the KPZ renormalization 
fixed point [20j|). Only recently we 21]. and other workers 



22H25|. were able to directly solve the continuum prob- 
lem and, until now, only within the droplet initial condi- 
tion. In addition to the convergence to F^is) at large t, 
the unveiled remarkable feature is that a proper gener- 
ating function g(s) (recalled below) remains a Fredholm 
determinant for all times t, hence leading to an exact so- 
lution for the universal crossover in time in the continuum 
KPZ equation. This universal distribution (depending on 
a single parameter t) describes in the DP framework the 
high temperature regime |2l| that has remarkable univer- 
sal features 26] . In the growth problem, this corresponds 
to a universal large diffusivity-weak noise limit, at fixed 
correlation length of the noise. 

In this Letter we obtain the corresponding exact result 
for the continuum KPZ equation for the case of flat initial 
conditions, most often encountered in experiments 0. 
We obtain here the generating function of the integer 
moments of the DP partition sum Z = e^ h (see below) 
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9^) = H-^r-^ > A = i(c 2 i/T 5 ) 1/3 , (2) 

n=0 



with c = DXq, T — 2v, as a Fredholm Pfaffian for any 
time t. The DP free energy, and the height field at a 
given point, take the form 
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h = hxZ = vot + \£ t 



(3) 



where g\(s) — exp(— e A ^*~ s )). In the large time limit 
goo(s) = limA^oo 9x( s ) = Prob(£ t < s) and we find that 
the distribution variable £t converge to the GOE Tracy 
Widom distribution Prob(£ t < s) = Fi(s). As in ^ 
we use the Bethe Ansatz (BA) for the replicated boson 
model with S attraction, and sum over all excited states, 
treating now the case of a DP with one end fixed and 
the other free. The calculation is far more complicated 



2 



than [2l| as we need the spatial integrals of the Bethe 
wavefunctions. Technically, this is surmounted by solving 
a half-space model in the proper limit. 

The solution of ([T]) for a given initial condition can be 

written, using the Cole-Hopf transformation Z : 



e 21/ n 



e 2 " 



£ fc 0M) = dyZ(x,t\y,0)e^ h( -y^ 



(4) 



in terms of the partition function of a DP, at temperature 
T = 1v and in the random potential V(x, t) = Xof](x,t), 
i.e the sum over paths x(r) G R starting at a;(0) = y and 
ending at x(t) = x 

Z(x, t\y, 0) = [ X(t) ~ X Dx(r)e-* X ^(£f + v( x (r),r)] ^ 

Jx{0)=y 

(5) 

with initial condition Z(x,t = 0\y, 0) = 5(x — y). In this 
continuum model the disorder correlation length is zero, 
i.e. V(x,t)V(x',t) = cS(t - t')5{x - x') with c = DXq, 
hence one can perform a rescaling x —> T 3 x/c and t — > 
2T 5 t/c 2 , and work in units such that T = 1 and c 
1 . as we do below. The continuum model (JSJ describes 
the high T limit of the DP on a lattice as discussed in 



21, 261. For the KPZ interface this continuum model 



describes the universal limit where the characteristic time 



t* = 2(2v) 5 /D 2 X% and space x* 



scales are much 



larger than the correlation lengths of the noise. 

In [2l[ we obtained the distribution of In Z for a DP 
with the two ends fixed Z = Z(0,*|0,0). From © this 
corresponds to a wedge initial condition for the KPZ in- 
terface, ^j/i W edge(£, t = 0) = — tu|a;|, in the limit of a 
narrow wedge w — > oo. Here we solve the flat inter- 
face initial conditions for KPZ, i.e. h(x,t = 0) = 0, 
hence the opposite limit w — > + of the wedge. It corre- 
sponds to the DP with one end fixed, one end free, i.e. 
Z&a,t{x,t) — J dyZ(x,t\y,0). To achieve that we found it 
easier to study the (left) half-space problem 



Z w (x,t) = / dye wy Z{x,t\y,0), 



(0) 



with Z(x,t = 0) = 6{-x)e wx . It is easy to show [31| 
that at w = + the half space model interpolates be- 
tween (i) the narrow wedge initial condition for x — > +oo, 
since the polymer is stretched and (ii) the flat (full-space) 
initial condition for x — > — oo, hence we study below 
limx-^-no lim,„_>.n Z,Jx, t) = Z Rat (x,t). 

As well known [27[ the calculation of the n-th integer 
moment of a DP partition sum can be expressed as a 
quantum mechanical problem for n particles described 
by the (attractive) Lieb-Liniger Hamiltonian [28[ 



H n = -1^-^- 25 1^ SiXi-Xj). (7) 
j — 1 3 l<i<j<n 

Generalizing (2lj . the quantum mechanical expectation 
for Z w (x,t) n is written as a sum over the un-normalized 



eigenfunctions \I/ Al (of norm denoted ||/x||) of H n with 
energies [33j : 

Z w (x,t) n = [ e w ^"=^( yi ...y n \e- tH -\x...x) 



-tE u 



W r 

*p := / e w ^y^^ yi ,..y n ), 
J Vi <o 



(8) 
(9) 



where we used the fact that only symmetric (i.e. bosonic) 
eigenstates contribute. The Bethe states ^ are super- 
positions of plane waves [28j over all permutations P of 
the rapidities Xj (j = 1, ..n) and we use the convention 



(10) 



where the coefficients A P = U n >i>k>i( 1 + - a^-a^*" )- 
The general eigenstates are built by partitioning the n 
particles into a set of n s bound-states formed by rtij > 1 
particles with n — X)j=i m j- Because we work with 
w = + , we can take directly the system size L = oo 
and in that limit (29| each bound state is a perfect string, 
i.e. a set of rapidities A J ' a = kj + ^(nij + 1 — 2a), where 
a = 1, mj labels the rapidities within the string. Such 
eigenstates have momentum = Y^j=i rn j^j an0 - en ~ 

ergy £; m = Y,%i{ m jtf ~ T5 m i( m i _ Tne ground- 
state corresponds to a single n-string with fcj = 0. 



In flHJ one already knows ^(x, .,x) - 
the norms of the string states [3C 

(c) r 



i!e-"^ A ° and 



A* = 



n\(Lc 



n-7 2 n 



3=1 l<i<j'<n 8 

(Jfei - fy) 2 + (m< + »"i) a c74 ' 1 J 



The new difficulty, i.e. computing the spatial integral 
of the Bethe states, simplifies dramatically for the half- 
space model. Using the symmetry of Eq. (fTU|). we have 
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Ai,..A„ 
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ju> + iAi + .. + i\j 



From the remarkable properties of the BA, it can be re- 
expressed, for any n and set of rapidities (with c = 1): 



ns=i(Ac 



n 



i + \ a + Xp — 2iw 



«W) X a + A/3 - 2iw 

l<a<p<n ^ 



(12) 
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If we now inject the string solution Aj 



Am, 



2a) + kj , we find after some elementary manipulations 



* lt =»!(-2rn^ ifc4 n 

i— 1 X<i<j<n s 

t(k u - m,-) 



D 



rrii.ki .rrij .kj 3 



r(Atii) 
r 



(13) 



D 



mi-\-m 3 



r 



with 



-iki — ikj — 2w + 1. 



We have now all ingredients to compute the generat- 
ing function <j2j) with Z = Z w (x,t). Writing the sum 
over states in © as all partitioning of n particles into 
n s strings and using that for L — > oo the string mo- 
menta rrijkj correspond to free particles [3(| (i-e. J2k ~~ ^ 
nfijL J = rrijL f k .), Eq. ([2]) becomes a sum over 



string configurations g\(s) = 1 + Z)^=i ^T-^( n s, s) with 



Z(n s ,x) = II 

m i , . . . m. n s —lj — l 



_2 m - 



nil 



n 

X<i<j<n s 



(14) 



Upon inspection we find that the limit of interest, 
w = + , is dominated by poles in the ki integra- 
tions, with S™. ik . ~ r(m ^k-+2w) and D %H,k i ,m j ,l h ~ 
i(k~+k )+2w ' ano - t na k the regular parts do not con- 

tribute. Replacing l/(ik + + ) — > ir5(k) yields an re- 
independent result, which can be shown to equal the limit 
x — > — oo, i.e. the flat initial condition for KPZ, on which 



we now focus. The result is the sum of the residues asso- 
ciated to configurations where the n s = 2N + M strings 
split into N pairs of strings of opposite momenta with 
same particle number m and M single strings of zero 
momentum with all distinct number of particles. After 
some non- trivial manipulations, detailed in [3lj ]. we bring 
the result in the form of a Pfaffian: 



%)=En| n 



^4j-m^ —Arrij kj A — Xrrij s 



xPf 



^fJiJk^Ji^ + q 

^S(h + k 5 ){-l) m ^ mitmj + i(2^) 2 ( 5(A; l ) ( 5(fc J )(-l) min(m - mj) sgnK - m s ) |(27r)<5(fc i ) 



(15) 



liki +mi -\-2ikj -j-rrij 



r 



We recall that for an antisymmetric matrix A of size 2n s arrive at our main result for g\(s) as a Fredholm Pfaffian 

(17) 



Pf^4= C- 1 )" II A CT (2i-i), CT (2i) , (16) 

<T£S 2ns ,cr(2j-l)<a(2j) i=l 



ffA ( s ) = Pf[J + K]= —,Z{n s ), 

n s — 

Z(n s ) = JJ / Pf[K(u i ,W J )] 2 „ s ,2n s 



with (PfA) 2 = detA We can now use the Airy trick 
[HEl e^ 3 ™* 3 = Uj I y MVj)e XV] and decouple the 
denominators in the lower right corner using auxiliary 
integrals fj ■ J v >Q e~ VjAi = pj ■ -j- and the numerators 

using derivatives. We perform rescaling kj — > kj/X, and . T / I \ . . , . , . . . . 

. ... .,9 m, , . ,i where J = T „ and K is an antisymmetric 2 by 2 

shifts yj —> yj + Vj — 4kj + s. Ihe summations over the \ —1 J 

rrii can be performed exactly inside the Pfaffian and we matrix kernel of components K a t = K a b(vi,Vj) with 
I 



K xx = I Ai( yi + v, + s + Ak 2 )Ai(y 2 + Vj + s + Ak' 2 )[ 



-2i(vi-Vj)k -rrJ\(b\ 

/ fc/A (e^+^)) + ^F(2e^,2e^)] 



2ik 
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Kl2 = \j M( - V 



Vi)(e~ 



1) 5(vj 



K 22 =28\v i -v j ), 

and the functions 

-2irkziF 2 (1;2 



h{z) 



2ik,2 + 2ik;-z) 



sinh (27rfc) T (2 
Zj) = sinh(z 2 



2ik)T(2 + 2ik) 



(18) 



(19) 



du 



x Jo(2a/z 1 z 2 (1 — u))[zi sinh^u) — z 2 sinh(z 2 'w)] . 

The full analysis of this result is performed in [3l|. Here, 
we first point out the simple one-string contribution 
{n s = 1), Z(l) = J v>Q K 12 {v,v) leading to 

Z{1) = \ [dy(e- 2 ° Xs -l)Ai(y + s), (20) 



also obtained [3l| from the ground state for each n, 
which gives the leading asymptotics of g\(s) for large 
A B 



s > 0. Using det 



C D 



= det I? det [A - BD- l C], 



g\(s) = Det [I — JK] can be written (for any time) in a 
form suitable for numerical evaluation [31]. The Pfaffian 
reported above allows simpler analytic manipulations. 

In the large time (large A) limit, one already sees from 
(pOf that Z(l) ->■ - / Ai(2y + s) = —TrB s where 
B s = 8(x)Ai(x + y + s)6(y) is the GOE kernel, as shown 
by Ferrari and Spohn [32j. This extends to all n s i.e. we 
find that 



lim Z{n s ) = (-1Y 



det[B s (xi,Xj)]„ a 



(21) 



Hence g-x>(s) — F 1 (s) — det[I — B s ] the Fredholm deter- 
minant expression for the GOE Tracy Widom distribu- 
tion. This is obtained using hniA-j.+oo fk/\i eXy ) = ~6{v) 
and limA-H-oo^e^, 2e^) = e(y 1 +y 2 ){8(y 1 )6(-y 2 )- 
Q{y2)9{—yi))- We checked (1211) e xplicitly up to n s = 4, 
but we report the proof in [3 if . For n s even it fol- 
lows from a slight generalization of [HI , namely det (I =F 
B s )/det(I±Bs) = / B>0 (J±B a )- 1 (ar s 0). 

To summarize we have obtained the generating func- 
tion for the distribution of the free energy of the DP with 
one free end, i.e. of the height of the continuum KPZ in- 
terface, for arbitrary time. At large time the distribution 
crosses over to the GOE Tracy Widom distribution F\ (s). 
Further properties of the finite time, including extracting 



P(f) and numerics are studied in [31 
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